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Two-step build-up of a thermoreversible polymer network: 
From early local to late collective dynamics 
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We probe the mechanisms at work in the build-up of thermoreversible gel networks, with the 
help of hybrid gelatin gels containing a controlled density of irreversible, covalent crosslinks (CL), 
which we quench below the physical gelation temperature. The detailed analysis of the dependence 
on covalent crosslink density of both the shear modulus and optical activity evolutions with time 
after quench enables us to identify two stages of the physical gelation process, separated by a tem¬ 
perature dependent crossover modulus: (i) an early nucleation regime during which rearrangements 
of the triple-helix CL play a negligible role, (ii) a late, logarithmic aging one, which is preserved, 
though slowed down, in the presence of irreversible CL. We show that aging is fully controlled by 
rearrangements and discuss the implication of our results in terms of the switch from an early, local 
dynamics to a late, cooperative long-range one. 

PACS numbers: 82.70-y,61.20.Lc, 61.43-j,62.20.F- 


I. INTRODUCTION 

Disordered solids, namely various glasses, polymeric 
or not, exhibit, following fast relaxation after quench, 
a slow, usually quasi-logarithmic, mechanical strength¬ 
ening. This self-decelarating behavior is known as aging. 
As pointed out by Parker and Normandjl), this specificity 
of glassy dynamics, commonly assigned to the complexity 
of the free energy landscape, is shared by another class 
of materials, namely thermoreversible hydrogels. Indeed, 
although their polymer content is in general very small 
(typically a few percent), when they are quenched be¬ 
low the gelation temperature T g , as crosslinks proliferate 
beyond percolation, the corresponding topological con¬ 
straints increasingly decimate the paths of easy relax¬ 
ation. 

While logarithmic aging has been evidenced in numer¬ 
ous gels 3, extensive systematic studies have focussed 
on gelatin ones. Gelatin consists of polypeptide chains 
which self-assemble, when quenched below T g ~ 30°C, 
via partial renaturation of the native collagen triple-helix 
structure m from the high temperature free chain coil 
configuration. The chirality of the resulting crosslinks 
(CL) permits to characterize this gelation dynamics by 
monitoring the evolution, not only of the storage shear 
modulus G, but also the optical activity, hence the rena- 
tured collagen fraction x [a]- 

There is general agreement n on the fact that gela¬ 
tion proceeds via the following scenario: 

(i) Formation of triple-helix nuclei whose critical size 
decreases as temperature T decreases below T g \ 

(ii) subsequent extension (growth) of these critical 
crosslinks while nucleation gradually slows down as the 
network mesh size decreases. 


The associated depletion of the entropy supply of the 
polymer chains leads to gradual exhaustion of these nu¬ 
cleation and growth processes, hence to a severe slowing 
down of the growth of G. This is often described as 
“saturation of the modulus”, since many authors tradi¬ 
tionally resort to the log G vs. t plot an example of which 
is shown on FigJT] (insert). However, in this late regime, 
as illustrated on Fig [1] 

(iii) structural relaxation does not stop but crosses over 
to the self-decelerating, logarithmic dynamics (G ~ logf) 
characteristic of glass aging. No trend toward saturation 
has been observed up to two months [7j ■ 



FIG. 1. Evolution of the shear modulus of a 5 wt.% gelatin 
gel in water at T = 20°C, vs. logarithm of time after quanch. 
The dashed line fits the late aging behavior. The double bar 
indicates the crossover region as deduces from optical data 
(see text section III.A.). Insert: same data in the lin-log 
representation. 


This aging behavior, specific of physical (as opposed to 
covalent) gels, has been assigned to network rearrange¬ 
ments taking place via thermally activated zippings and 
* tnstan@msp.jussieu.fr unzippings of the H-bonds stabilized triple helices, i.e. to 
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the thermoreversibility of the CL. However, no direct in¬ 
formation about the relative weight of nucleation/growth 
versus “crosslink plasticity” in the early and late regimes 
of network build-up has been available up to now. In 
order to try and shed more light on this issue, we resort 
in this work to an approach which parallels the study 
of point-to-set correlations in glass-forming liquids |8j. 
That is, in order to assess the influence of reversibility 
on the dynamics, we introduce in the gels a controlled 
density of covalent CL. 

More precisely, we first form at a temperature above 
T g a covalent network, characterized by its shear modu¬ 
lus G c h■ Further covalent crosslinking is then inhibited, 
after which we quench the gel below T g so as to allow 
for physical crosslinking to take place. This leads to the 
formation of “chemical-physical hybrid gels” @, HI the 
aging dynamics of which we monitor via the evolution of 
their modulus G(t) and renatured collagen fraction \(t). 

As can be intuitively expected, increasing G c h , hence 
the fraction of irreversible CL, results in a slower and 
slower dynamics. Moreover, the qualitative features 
of x{t) and G(t) — including non-saturating, logarith¬ 
mic aging of the shear modulus — are fully preserved. 
Detailed quantitative analysis of the dependence upon 
G c h of network build-up dynamics enables us to iden¬ 
tify a temperature-dependent cross-over shear modulus 
separating (i) an early nucleation regime during which 
crosslink rearrangements play a negligible role, from (ii) 
the late, logarithmic aging one, fully controlled by rear¬ 
rangements. We discuss the implication of our results in 
terms of the switch from an early, local dynamics to a 
late, cooperative long-range one. 


II. MATERIALS AND METHODS 

Sample preparation — The first stage consists in the 
preparation of a covalently bonded gelatin network with 
a given modulus G c h■ For this purpose we use an enzy¬ 
matic route which we have described in detail elsewhere 

Q 

In short, we dissolve gelatin (300 Bloom, type A from 
porcine skin, Sigma) in deionized water et 65° C. After 
total dissolution of the polymer, the solution is quickly 
mixed at 40°C with a Tgase enzyme solution (microbial 
transglutaminase, Activa-WM, Ajinomoto Foods Europe 
SAS) so as to reach a final composition of 5 wt% gelatin 
and 2.6 nmol of Tgase (corresponding to an enzymatic 
activity of 2U). The solution is poured into the mea¬ 
surement cells of both a rheometer and a polarimeter 
(see below). Gelation then proceeds at T set = 40°C. 
At this temperature, chosen well above T g cs 30° C, no 
triple-helix reversible cross-link can form whereas Tgase 
catalyzes actively the formation of inter-chain covalent 
bonds between two specific residues [l2|. 

When the target shear modulus level, as monitored 
in the rheometer, is reached, the enzyme is inhibited 
by quickly heating both samples up to 70°C and stay¬ 


ing at this temperature for 10 min, after which they are 
cooled down back to 40°C. We have checked that the 
shear modulus, measured at this temperature, remains 
constant over days. 

The physical gelation is subsequently triggered by 
quenching to the working temperature T <T g in most 
cases, T = 20°C. 

Shear modulus measurements — The elastic response 
of the hybrid gel is characterized in a stress-controlled 
rheometer (AMCR 501, Anton Paar) equiped with a 
plane-plane, sand-blasted cell, oscillating at 1 Hz, with 
a 1% strain amplitude so as to operate in the linear re¬ 
sponse regime. The sample is protected against solvent 
evaporation by a paraffin oil rim. We quench the sam¬ 
ples at the maximum rate (15°C/min) permitted by the 
thermoelectric cooling device so as to minimize the for¬ 
mation of physical crosslinks prior to reaching the target 
temperature T. 

We will need to characterize the hybrid gels by the 
modulus of the purely covalent network at temperature 
T. Since physical gelation starts without any time-lag 
(see below), it is not possible to measure reliably G c h{T). 
Since the gel elasticity is of entropic origin, the shear 
modulus at fixed network structure scales as the absolute 
temperature T. So, as a proxy for G c /j(T), we use the 
following expression: G ch {T) = (T/T set )G ch (T set ) with 
both temperatures in °K. 

All the gels studied here exhibit a very small level of 
viscoelastic losses (typically the ratio of the loss to stor¬ 
age moduli, G"/G' < 10 -2 at 1 Hz). So, for the sake 
of simplicity we characterize them by their sole storage 
shear modulus, noted G. 

Optical activity measurements — The rotation angle 
a of the polarization plane of light passing through a 1 
cm gelatin sample is measured using a spectropolarime- 
ter (Jasco 1100) operating at 436 nm. Gelatin is optically 
active in both the random coil (above T g ) and helix (be¬ 
low T g ) conformations. The renatured collagen fraction 
X is obtained as [1, [§]: 


/ '- Q,collagen _ Q,coil ' ' 

where a C011 is the rotation angle measured in the coil 
state above T g and extrapolated down to T according to 
Q, and a colla g en is the rotation angle of a 1 cm-thick, 5 
wt.% collagen sample 

The sample cell is thermalized by a water circulating 
circuit. Quenching is achieved by switching the circulat¬ 
ing water between two different temperature-controlled 
baths. Although the total duration of the quench is com¬ 
parable to the one achieved in the rheometer, it is impos¬ 
sible to match precisely the shapes of the transients, the 
initial cooling rate being faster in the polarimeter. In 
the following we will chose the time origin for physical 
gelation as the instant at which the sample temperature 
reaches the target one T. 
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III. RESULTS AND DISCUSSION 




FIG. 2. (a) Shear modulus dependence on time after quench 
to T = 20° C of a set of hybrid gels with various densities 
of covalent crosslinks. The labels indicate the values of the 
shear modulus G c h of the corresponding chemical networks, 
(b) Evolution of the renatured collagen fraction \ for the same 
gels. Data from the G c h = 1480 Pa gel have been discarded 
due to the smallness of the corresponding signal to noise ratio. 


Figure [2] shows the time evolution of the shear modulus 
G and helix fraction x for values of the modulus of the 
covalent network at the working temperature T = 20°C, 
increasing from 0 (purely physical gel) to 1480 Pa. The 
semi-logarithmic representation (FigJSJa) reveals that the 
evolution of G is qualitatively preserved in the presence 
of covalent CL. In particular, late quasi-logarithmic aging 
is observed for all values of G c h , no trend to saturation 
being observable up to 10 5 s. 

The early evolution of the modulus is best displayed by 
plotting the increment AG of its value from that immedi¬ 
ately after quench. The most conspicuous feature of the 
corresponding plot (FigO consists of the initial behav¬ 
ior of AG: while the purely physical gel exhibits a finite 
delay before the emergence of a measurable modulus, for 
all hybrid gels, however small G c h > 0, the initial rate of 
growth of AG is finite, decreasing as G c h increases. The 
origin of this discontinuous behavior is clear: while, in the 
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FIG. 3. Early evolution of the excess shear modulus A G(t) = 
G(t) — Gch■ Same data as Fig[2]a. 


purely physical gel case, percolation of the triple-helix CL 
has to be reached for a finite G to emerge, in the pres¬ 
ence of a preexisting percolating covalent network, phys¬ 
ical CL formation is immediately revealed. This leads us 
to an important inference, namely a fraction at least of 
the triple-helix CL involve polymer strands belonging to 
the covalent network. This proves that our hybrid gels 
do not consist (as is the case of “double network” ones 
[H) °f t wo interpenetrating but independent networks of 
different natures — which would contribute additively to 
G. Finally, the decrease of the initial slope indicates that 
covalent interchain bonds do not act as inhomogeneous 
nucleation centers for triple helix formation. 

The evolution of the helix fraction \ is shown on 
FigOJb. Here again the hindering effect of covalent CL 
translates into the gradual slowing-down of the dynamics 
with increasing CL density. Unlike G, the optical activ¬ 
ity of the purely physical gel does not exhibit any latency 
since it is insensitive to the connectivity of the polymer 
chain system. 

Beyond these general remarks we need to define a 
procedure to compare modulus build-up in our various 
gels. We choose to proceed as follows: we first measure 
the time to(G c h) for which the modulus of the purely 
physical gel G p h{t) reaches the value G c ^, then we shift 
the G(t; G c h) curve of the corresponding hybrid gel by 
to{G c h ) along the time axis. 

The rationale for this choice is the following. A purely 
physical gel and a chemical one exhibit distinct network 
topologies, since triple-helix CL are spatially extended, 
involving three chains, while the pointlike covalent ones 
connect two chains only. However, as gel elasticity is of 
purely entropic origin, G is controlled by the entropy re¬ 
duction of the polymer chain system due to crosslinking, 
hence is an indicator of the density of degrees of freedom 
which are frozen at the corresponding stage of gelation 
Q So, our procedure amounts to choosing for the com¬ 
mon reference the situation where the density of degrees 
of freedom which remain available for further evolution 
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for both gels is comparable. 



FIG. 4. Same data as Fig[2]a .: each hybrid gel curve has been 
shifted along the time axis so as to allow comparison with the 
purely physical gel one G p h{t) (black curve, full dots). to{G c h) 
is the time at which G p h = G c h . The hatched strip indicates 
the end of the early regime. The double vertical bar indicates 
the cross-over range as obtained from x(G) data (see text 
section III.A.) 

The resulting plot (Fig|4j exhibits a striking structure. 
Namely, up to a modulus level G* lying in the 700-800 
Pa range, the hybrid gel data corresponding to G c h < G* 
closely bunch with G p h{t). Beyond this, they separate 
from Gph and bend sharply into a splayed set of quasi 
log-linear curves. For G c h > G*, separation occurs for 
G — G c h values which decreases rapidly with G c h ■ 

This structure leads us to identify, in the network 
build-up process, two clearly different regimes separated 
by a cross-over: 

(i) an early stage in which the dynamics is basically 
insensitive to the nature of preexisting CL. 

(ii) a late, aging one in which the modulus growth is 
strongly affected by the degree of CL reversibility as re¬ 
vealed by the marked decrease with increasing G c h of the 
logarithmic slope. 

We now discuss separately in more detail the features 
of these two regimes. 

A. Early stage 

The bunching of the G curves (FigQJ suggests that, in 
this regime, it is the total gel modulus which codes the 
network evolution. If such is the case, other characteris¬ 
tics of the network should exhibit a change of regime at 
a the crossover value comparable to G*. A natural can¬ 
didate for checking this proposition is the helix fraction 
X which we therefore plot, in the “master curve” repre¬ 
sentation of Joly-Duhamel et al. [If|, as a function of 
G. 

FigEJa shows the corresponding plot for a purely phys¬ 
ical gel at T = 20°C. The finite x(G = 0) value measures 



FIG. 5. (a) Helix fraction \ vs. shear modulus G for the 
physical gel at T = 20°C. The thick dashed line is the fit to 
the initial linear regime ending at G = G* x . Insert: Inverse of 
the linear regime slope, dG/d\, vs. reduced departure from 
gel temperature T g . (b) Helix fraction versus shear modulus 
for a hybrid gel with G c h = 253 Pa at T = 20°C. Insert: End 
of linear regime G x vs. chemical modulus G c h at the same 
temperature. 


the helix fraction at percolation. At the modulus value 
noted G* = 770 ± 40 Pa, x(G) crosses over from a linear 
regime to a markedly sublinear one. The same analysis of 
the hybrid gel data for G c h values up to 750 Pa reveals, as 
illustrated on FigEJb, the persistence of a low-G linear 
regime terminating at G*(G c h). We find (see insert of 
Figj5]b) that G* does not exhibit any significant depen¬ 
dence on G c h■ The mean (± std. dev.): G* = 830 (±90) 
Pa is fully compatible with the crossover level G* previ¬ 
ously identified from the bunching of the G(t) data (see 
Fig|4j. Moreover, the G* crossover turns out (Fig{T]) to 
lie very close to the emergence of the logarithmic aging 
regime. In order to check the robustness of these fea¬ 
tures with respect to temperature we have performed the 
same analysis on physical gels of the same composition, 
quenched at T = 15 and 10°C. As is well known, the 
smaller T, the faster the gelation dynamics. The shape 
of the x(G) curves is preserved while G* increases as T 
decreases — its values, reported in Table U , remaining 
in the close vicinity of the entry into the log stage (see 
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Fig®. 


from 20 to 1600 Pa. These data are very nicely fitted by 
the expression 


TABLE I. default 


T [°C] 

<5; [Pa] 

G [Pa] 

10 

2650 

2016A310 

15 

1670 

1300A160 

20 

830 

685±15 


These results bring strong support to our above as¬ 
sumption of the existence of an early regime, the devel¬ 
opment and termination of which are coded by the to¬ 
tal elastic modulus, i.e. are quasi-insensitive to network 
topology and crosslink reversibility. 

Moreover, the quasi-linearity of the early \{G) de¬ 
pendence suggests a simple, qualitative interpretation, 
namely: according to the entropic origin of gel elasticity, 
it is commonly admitted fltl that, as a first approxima¬ 
tion, the shear modulus scales as the CL density. During 
physical gelation, only triple helix CL are formed. So, 
on the same level of approximation, we may write for 
the variations of G and \ induced by an increment Av 
of CL density: AG ~ ksTAv and Ay ~ tAv with t 
the average triple-helix CL length. From this, d\/dG = 
cst immediately translates into the statement that l re¬ 
mains stationary over the early gelation regime. Note 
that this assertion agrees with the conclusion reached by 
Guo et al 0 who were able to deduce, from the evolu¬ 
tion of x upon remelting physical gels of various ages, the 
distribution of lengths of these rods at various stages of 
gelation. They found that it remains quasi invariant in 
the corresponding time regime (while £ slowly drifts, at 
later times, toward larger values). 

One step further, in agreement with nucleation theory 
eh, we can expect the critical nucleus size, hence i to 
diverge upon approaching the gelation temperature T g . 
Since, from our above remarks, l/£ ~ dG/d\, we have 
plotted dG/dy vs. (T g — T)/T g . We find that our data 
for T = 10,15, 20°C extrapolate linearly to the origin 
(see insert of Fig® a), indicating the expected (T g —T) _1 
divergence. 


B. Late stage 


r = Fq exp 



( 3 ) 


with Q = 730 ± 15 Pa. Note that the extrapolated slope 
r 0 = 466 ± 4 Pa/decade lies noticeably below the value 
(630 Pa, see black dot on Fig® for the purely physical 
gel. This is attributable to the fact that, while the physi¬ 
cal gel does not contain any irreversible CL, a hybrid gel 
with vanishingly small G c h presents the finite density of 
covalent CL corresponding to the percolation threshold 
of the chemical network. The above exponential decay is 
robust with respect to temperature variations as shown 
on the insert of Fig® which reports data obtained at 
T = 10 and 15°C. 



G 

ch 


[Pa] 


FIG. 6. Logathmic slope F = dG/d log 10 t vs. G c h for hybrid 
gels at T = 20°C. Full line: best fit T = Aexp(—G c h/G) with 
G = 685 Pa. The purely physical gel data point (full dot) has 
been excluded from the fit (see text). Insert: T vs. G c h at 
different temperatures below T g . 


This behavior brings to light the existence of a charac¬ 
teristic modulus G(T) associated with logarithmic aging. 
We notice that its variation with T parallels that of the 
crossover modulus while its values remain systematically 
smaller than G* (see Table Hand Fig. Q. 


This regime is characterized by the logarithmic in¬ 
crease with time of the shear modulus which emerges 
close above the crossover level G* defined on Fig® As 
already mentioned, the logarithmic slope: 


d(logio t) 

is markedly dependent upon G c h, i.e. on the density of 
irreversible CL. Fig® displays the values of T measured 
at T = 20° C for various hybrid gels with G c h ranging 


C. Discussion 

These results demonstrate the existence of two dynam¬ 
ical regimes during the build-up of the physical network, 
which can be clearly differentiated according to their sen¬ 
sitivity to the nature of crosslinks. 

While covalent CL are irreversible, physical ones un¬ 
dergo “rearrangements” via elementary, thermally acti¬ 
vated zipping or unzipping events through which helix 
segments may either grow, shrink or slide by coherent zip¬ 
ping/unzipping at both ends. The clock is set for these 
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FIG. 7. Modulus build-up curves G p h{t) at different temper¬ 
atures. The circles indicate the cross-over modulus G£. The 
squares indicate the modulus Q characteristic of the sensitiv¬ 
ity of gel aging to G c h (see Tabic [I]) 


fluctuations by the cis-trans isomerisation of (Pro-Gly- 
Pro) sequences of amino acid residues d MM- As the 
corresponding mean flip time lies in the 10-100 s range 
at T — T gi we are dealing here with very slow processes. 
Rearrangements occur via the biasing of fluctuations un¬ 
der the effect of tension forces along the coil segments 
which connect the CL. 

In the early regime, the evolutions of both the stiff¬ 
ness and the helix fraction are found to be coded by the 
instantaneous value of the modulus and essentially insen¬ 
sitive to the presence of covalent bonds. This implies that 
the time scale for the nucleation of triple helices is small 
enough for rearrangements to be quasi-frozen, hence the 
insensitivity to the (reversible or not) nature of the CL. 

Of course, the larger G , the shorter the mesh size hence 
the thinner the entropic supply available for further nu¬ 
cleation the rate of which thus decreases — so much so 
that rearrangement dynamics finally takes over. This 
results in a rather sharp switch, about a crossover mod¬ 
ulus G*. The emerging late regime is characterized by 
the logarithmic growth of G(t) together with a markedly 
sublinear x(G) — which indicates that less and less new 
helix length is required for further stiffening of the net¬ 
work. Moreover, and most important, the aging dynam¬ 
ics dramatically slows down as the density of irreversible 
CL increases. All these features are consistent with a 
dynamics controlled by CL rearrangements. 

Loosely speaking we may distinguish between two 
types of such processes, which we refer to as local and 
collective. We define as a local process the length evolu¬ 
tion of a triple helix segment in the frozen environment 
determined by its non-evolving nearest neighboring CL. 
Due to the finite, nanometric, size of the gel mesh, the 
relaxation time for such a process is clearly finite. A 
simple model of the helix-coil transition on a single poly¬ 
mer strand with fixed ends !2fl] and with a Fokker-Planck 
dynamics [2l[ for the helix length evolution yields a lin¬ 


ear dependence of the relaxation time on the number of 
monomers in the strand f22| . Explaining the late dy¬ 
namics by the “local” growth of triple helices nucleated 
in the early stage would demand a flat distribution of 
relaxation times spanning more than the 3 decades over 
which logarithmic aging has been observed, hence a com¬ 
pletely unrealistic flat (gaussian) strand end-to-end dis¬ 
tance distribution spanning at least 1.5 decades. 

We are therefore led to conclude that the aging 
behavior results from dynamical interactions between 
crosslinks, which we call collective processes. More pre¬ 
cisely: an elementary (un)zipping event affecting a phys¬ 
ical CL induces variations of the tensions of emerging 
strands, i.e. acts as a localized force dipole acting on 
the surrounding elastic medium. This gives rise quasi- 
instantaneously (i.e. with a delay set by sound propaga¬ 
tion) to a power law decreasing strain field which, in turn, 
modifies the tensions of strands involved in remote CLs, 
thereby displacing their equilibrium triple helix length. 
Due to the long range of elastic couplings, this results 
in a broadband mechanical noise. As global relaxation, 
i.e. collagen renaturation, proceeds, one expects that 
further modulus growth requires cooperation between an 
increasingly large number of physical CL, hence involves 
a growing correlated volume and, accordingly, a growing 
time scale. 

In such a picture, introducing an increasing fraction of 
irreversible CL has a multiple effect. On the one hand, 
for a given G it reduces the fraction of rearrangeable CLs. 
On the other hand, the average distance between the lat¬ 
ter increases, and so does the mechanical noise ampli¬ 
tude. Moreover, covalent bonds restrict the amplitude of 
triple helix sliding due to coherent zipping/unzipping at 
both ends. We saw that as time lapses in the late regime, 
modulus increase involves less and less helix length cre¬ 
ation. This is a strong hint of the increasing role played 
by length preserving correlated sliding processes, which 
are destroyed by the interposition of covalent bonds on 
chains connecting the two helix segments. This topologi¬ 
cal effect of chemical bonds thus results in the decimation 
of the possible relaxation paths in configuration space. 

At this stage, pending the development of a theory 
relating logarithmic aging to definite characteristics of 
structural evolution, it remains out of the scope of the 
present work to explain on the basis of these qualitative 
remarks the observed dependence of the aging dynamics 
of hybrid gels on their chemical modulus G c h■ Progress 
toward this goal would certainly highly benefit from the 
development of numerical modeling, possibly in the spirit 
of the investigations of long time scale glass aging cur¬ 
rently developed by several groups [2314251 1 ■ 

This might shed light on an issue of broader interest, 
namely whether or not stiffening during aging is asso¬ 
ciated with decreasing elastic non-affinity, i.e. whether 
slow relaxation drives the system toward gradual homog¬ 
enization of the internal stresses 
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